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The effects of blending on the light curve shape of Cepheids 
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Abstract. A short analysis is presented of the effects on the cepheid light curve shape, i.e. on the Fourier parameters 
usually adopted for its description, of the blending of the stellar image with other close stars. The conclusion is 
that, within reasonable error, the effects are in general small and the Fourier decomposition is confirmed to be a 
useful tool for pulsation mode discrimination. A large effect has been found on the phase differences in a narrow 
period range corresponding to the known resonance centers between pulsation modes. 
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1. Introduction 

Cepheids are primary distance indicators for external 
galaxies and those used for this application pulsate in 
the fundamental mode. First overtone mode Cepheids are 
brighter by about 0.4 mag than fundamental mode pul- 
sators with the same period. Since the period-luminosity 
relation has an intrinsic dispersion, which depends on sev- 
eral parameters (e.g. different effective temperature or 
color, different reddening, contribution from stellar com- 
panions), it is essential to remove the contaminating stars 
that are pulsating in a different mode. The large surveys of 
the Magella nic Clouds performed by M ACHO (e.g. Welch 
et al. [L997D , EROS (e.g. Beaulieu et al. |1995| ) and OGLE 
(e.g. Udalski et al. 199S ) projects proved that the Fourier 
decomposition is a good technique for discriminating the 
mode among short period (P 6 d) Cepheids. More re- 
cently, the technique began to be applied to Cepheids 
of farther galaxies in the Local Group, such as IC 1613 
(e.g. Antonello et al. |1999| Dolphin et al. |2001[ ) and M33 
(Mochejska et al. pOOip . 

The large surveys offered also the opportunity of dis- 
cussing the problems related to blending. Mochejska et al. 
( 200C ) define the blending as the close projected associ- 
ation of a Cepheid with one or more intrinsically lumi- 
nous stars, which cannot be detected within the observed 
point-spread function by the photometric analysis. There 
is some debate about the implications for the distance de- 
termination related to the blending and more generally to 
poor resolution of the stellar images in these galaxies. The 
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blending also has other effects on the light and the color 
curves. Mochejska et al. ( |2000 ) note that in the case of 
a red or blue companion the light curve exhibits a flatter 
minimum. As regards binaries, it is well-known that the 
observed amplitude of the light curve is affected by the 
luminosity of a bright companion. Could it be that the 
blending, apart from producing a lower amplitude, also 
mimics a different pulsation mode? Recently, we recalled 
that in principle such an effect on the Fourier parameters 
is small in the context of mode identification (Antonello 
et al. 2002 ). Here we report the results of simulations that 
support this conclusion, and we discuss some unexpected 
characteristics. 



2. Analysis 

The problem by itself would not be important if we adopt 
intensities instead of magnitudes to measure stellar bright- 
ness. Indeed, an increased intensity due to a close star, as- 
suming no measurement error, would produce a light curve 
with a similar shape to that without such a close star. The 
average intensity would be larger, the absolute amplitude 
would be the same, and the relative amplitude would be 
of course decreased. Let < X > be the average stellar in- 
tensity (that is, the average number of collected photons), 
AX the absolute amplitude, A = AL/ < L > the rela- 
tive amplitude, and e ~ \J < L > the mean absolute error 
on the measurement. Let us assume a close constant star 
with intensity a < L >. The relative amplitude of the 
system will be A\ — AX/ [(a + 1) < X >] and the mean 
absolute error t\ ~ y/(a + 1 ) < X >. A close star has the 
effect of decreasing the relative amplitude and increasing 
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Fig. 1. Lower panel: blending effect on the V light curve 
of a Cephcid, for different values of the luminosity of the 
companion star (L c ). Upper panel: comparison between 
the light curve for L c = (continuous line) and L c = 4 < 
L > (dotted line) scaled to the same amplitude. 



the absolute error. This implies a lower order of fit of the 
reliable Fourier decomposition of the intensity curve, and 
larger formal errors of the Fourier parameters; however, 
the parameters themselves are unchanged (within the for- 
mal errors). 

The nonlincarity of the relation between intensity and 
magnitude introduces some changes. The simplest method 
for studying them is by means of simulations. We consid- 
ered light curves of some stars pulsating in the fundamen- 
tal or first overtone mode ( e.g. X Cyg, DT Cyg) observed 
by Moffett & Barnes ( 1984 ; data retrieved from McMaster 
Cepheid Photometry and Radial Velocity Data Archive), 
and we adopted the best fitting curve as a synthetic light 
curve. We simulated several time series, adopting the origi- 
nal observing dates, and changing the synthetic light curve 
by introducing the contribution of a close constant star, 
and different mean errors of the measurement. In Fig. [j]we 
show the effects of increasing luminosity on the synthetic 
light curve of X Cyg. In the upper panel one can see the 
changes of light curve shape due to a four times brighter 
companion; the two curves are scaled to the same ampli- 
tude. The flattening of the minimum does not appear very 
prominent, even in this case where the magnitude differ- 



ence between the Cepheid and the blended image is large, 
1.75 mag. 

The time series were constructed applying a random 
number generator for a Gaussian error distribution. The 
series were then Fourier decomposed and the resulting 
Fourier parameters are plotted in Fig. |^ for the case of 
X Cyg, as an example. One can see clearly that the in- 
creasing blending implies a decreasing order of the reliable 
fit. 

When performing the simulations, we also analyzed 
some OGLE stars in the SMC, and we noted different 
trends with respect to the above Cepheids. We suspected 
some dependence on the P, therefore we decided to an- 
alyze all the Cephe ids in OGLE database of the SMC 
(Udalski et al. |l99S|) . The fitting curves of the Fourier de- 
composed /-band light curves were modified by introduc- 
ing the contribution of a companion star with L c = 2 < 
L >, then they were analyzed and we computed the dif- 
ference between the Fourier parameters for L c = 2 < L > 
and L c = 0. The results for the lowest order are shown in 
Fig. H and ^ for the fundamental and first overtone mode, 
respectively. Although the effect on the amplitude ratio 
is always small, the trend with P is confirmed. The un- 
expected result is the large effect on the phase difference 
very close to the resonance centers at P ~ 10 d for the 
fundamental mode, and P ~ 2.2 d for the first overtone 
mode. Outside these narrow P ranges the effect is small. 

3. Discussion and conclusion 

The cases discussed here concern reasonable light curves; 
we do not consider the problems related to very faint vari- 
ables, which can hardly be detected at minimum light. The 
requirement is that in the P interval where it is possible 
to find stars pulsating in different modes, the Fourier pa- 
rameters must allow us to make the discrimination. It is 
known that this occurs for P ^ 6 d for the fundamen- 
tal and the first overtone mode, using only light curves 
parameters. The results of the simulations show that in 
this P range the blending has a negligible effect when we 
compare the differences introduced by it with the size of 
the parameters themselves. In particular, a blended funda- 
mental mode pulsator will have slightly larger amplitude 
ratios than a non-blended one; we recall that the ampli- 
tude ratios of fundamental mode pulsators are intrinsi- 
cally larger than those of first overtone mode ones in this 
P range. The same occurs for a first overtone mode pul- 
sator compared with a second overtone one, for P ^ 1.3 d. 
On the other hand, a heavily blended first overtone pul- 
sator increases its R21 value, but in general not so much 
so as to be confused with a fundamental mode pulsator. 
In conclusion, the blending due to various reasons is not 
an issue for the pulsation mode discrimination. 

The color of the companion stars is not relevant for 
the present discussion, as long as their contribution is 
constant; some (second order) effects could be related to 
their intrinsic variability, both in terms of photometric 
variability and/or Doppler shift. The influence of the pho- 
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Fig. 2. The plots show how the Fourier parameters and 
light curve amplitude of a Cepheid change according to 
the luminosity of a companion star [L c / < L > is the ra- 
tio of the luminosity of the companion to the average value 
of the Cepheid). The symbols indicate different values of 
the mean error a of measurements adopted in the simula- 
tions: filled circle: cr=0, filled triangle: cr=0.02; open circle: 
cr=0.05; open triangle: (7=0.1 mag. The errorbar indicates 
the formal error of the respective parameter. Amag is the 
average magnitude difference between the Cepheid and 
the blended image 



tometric variability of the companion itself can be usu- 
ally accurately estimated, since an adequate time series 
analysis is sufficient to disentangle the different contribu- 
tions, because of the different periodicities or timescales 
involved. Also in this case, however, it is wise to work with 
intensities rather than with magnitudes. Variable seeing 
conditions could have some effect on the estimate of the 
intensity through the PSF fitting procedure; however in 
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Fig. 3. Simulated blending effect on the /-band light 
curves of all the OGLE fundamental mode Cepheids in 
the SMC. The plots show the difference of R21 and fai 
between the light curves for L c = 2 < L > and L c = 



this case we would expect just an increased error in the 
measurement. 

The plots in Figs. || and |J suggest some interesting 
considerations. A light curve with an altered value of the 
mean luminosity, such as that depicted in Fig. 1, or ex- 
pressed with a different, nonlinear mathematical function 
(e.g. the intensity instead of the magnitude) is character- 
ized of course by (usually slightly) different Fourier pa- 
rameters. If we estimate the differences related to these 
changes, we note that the largest ones are for the phases 
of the Fourier components with smaller amplitude; for ex- 
ample, at about 10 d some stars have R21 < Rn, for i 
from 3 up to 6 or more. The large differences are not due 
to errors or to uncertainties, since here we are not dealing 
with observed data but with synthetic light curves (i.e. the 
fitting curves), which are in principle error-free. In other 
words the differences are intrinsically real and reflect di- 
rectly the change of the shape introduced by the different 
mathematical function. The interpretation of this feature 
is reported in the Appendix; from that, we conclude that 
the observed dispersion is strictly related to the smallness 
of the Fourier component involved. In our example, the 
small second Fourier component has changed its phase 
value by several tenths of a radian, while for the other 
components the change is much smaller. For the same rea- 
son we should expect an analogous results for ^41, i.e. we 
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Fig. 4. Simulated blending effect on the J-band light 
curves of all the OGLE first overtone mode Cepheids in 
the SMC. The plots show the difference of Ri\ and </>2i 
between the light curves for L c = 2 < L > and L c = 



Fig. 5. The difference of R41 and (f>4i between the sim- 
ulated light curves for L c = 2 < L > and L c = of 
fundamental mode Cepheids. It should be compared with 
Fig. 1 



should have some dispersion at P ~ 7 d, where R41 is 
small since another resonanc e, P0/ P4 = 3, should be op- 
erating there (e.g. Antonello 1994 ). Indeed this is shown 
in Fig. ||; note also that the discontinuity of Ai? 2 i located 
at 10 d is replaced by that of Ai?4i at about 7 d. In a 
certain sense, plots such as those shown in Figs. [|, || and 
H are better indicators of resonance effects than the classi- 
cal ones, because they are free of subjective corrections of 
the phase differences by ±2n, which could be uncertain, 
mainly for the higher orders. Finally, it is possible to note 
two minima in the lower panel of Fig. [| one at the reso- 
nance center, and the other at logP ~ 1.5. Kanbur et al. 
(2002) noted the structural change of the light curves at 
this P; these features still await a theoretical interpreta- 
tion. 

Last but not least, we remark further that several prob- 
lems with the time series analysis of stellar luminosities 
would be simplified by adopting intensity scales instead 
of magnitude scales. This statement is not new, of course. 
Our comment is just further support to the proposal of 
abandoning the magnitudes. In fact, the blending has no 
effect on the light curve shape when we use intensity light 
curves, and this is an advantage, since one is always deal- 
ing with observed parameters which are affected by errors. 
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Appendix A: Intensity and magnitudes 

Note that the increasing blending produces a light curve, 
expressed in magnitudes, with a shape which is similar 
to the shape of the intensity light curve. That is, for 
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very large L c , the amplitude becomes very small, and 
the Fourier parameters become those of the intensity- 
light curve (for the phase differences one has to consider 
the different sign). The diagrams of the differences be- 
tween intensity- and magnitude-light curves of the SMC 
Cepheids look similar to those of the diagrams shown in 
Figs. 3, 4 and 5, but with slightly different ranges of the 
ordinatae; for example, the range of Ai?2i values would 
be about ±0.05 instead of about ±0.035 as indicated by 
Fig. 3. 

In this Appendix we will use some approximations to 
understand the effect seen near the resonances between 
pulsation modes of Cepheids, or more generally the ef- 
fect on the smaller Fourier components, given by different 
mathematical descriptions of the light curve. In this con- 
text, the intensity-light curve could be considered, in a 
certain sense, as a magnitude-light curve for an extremely 
large blending value. 

Let us assume that the intensity curve is expressed by 

L = L Q + x = L a + J^[^4i cos(iojt) + Bi sin(iwt)], (A.I) 

where L Q is the mean intensity value, which may include 
the contribution from a companion star or blending, and 
to is the pulsation frequency. The observed light curve can 
be written as 

V = -2.51og(L) + fci, (A.2) 

where k\ is an appropriate constant. By considering the 
natural logarithm, we can write 

V = ln(L) + k[ = ln(l + x/L Q ) + k 2 , (A.3) 

where V = -V/1.0857 and k 2 = -fc 1 /1.0857+rn(£ o ). We 
assume a relatively small amplitude, and expand (A.3) in 
the series 

V = k 2 +x/L -(x/L ) 2 /2 + ... (A.4) 

where the Fourier series of espression (A.l) is introduced, 
and we assume for simplicity that i < 3. After some ma- 
nipulation, we get the following expressions for the cocffi- 



and the correcting term for the mean value: 
-{A\ 



B\ 



Bi + Ai + Bi)/ALi. 



(A.17) 



cients of the cosine terms, from i 


= 1 to i = 6, 




At/L - {A X A 2 + B 1 B 2 + A 2 A 3 - 


f B 2 B 3 )/2L 2 Q 


(A.5) 


A 2 /L -(A 1 A 3 +B l B 3 +A 2 /2 


+ Bf/2)/2L 2 


(A.6) 


A 3 /L a - (A X A 2 - 


- B X B 2 )/2L 2 


(A.7) 


-{A X A 3 -B X B 3 + A\j2 


-Bl/2)/2Ll 


(A.8) 


-(A 2 A 3 - 


- B 2 B 3 )/2L 2 


(A.9) 


-(At/2- 


- B 2 3 /2)/2L 2 , 


(A.10) 



respectively. Six Fourier components are needed instead 
of just three to describe the V light curve. Analogously 
for the sine terms we get 

BxlL - (A X B 2 - B X A 2 + A 2 B 3 - B 2 A 3 )/2L 2 Q 

B 2 /L - (Ai-Bi + A X B 3 

B 3 /L - (A X B 2 

-(Arfs + BtAs 

~(A 2 B 3 



- B X A 3 )I2L 2 
+ B 1 A 2 )/2L 2 a 
+ A 2 B 2 )/2L 2 a 
+ B 2 A 3 )/2L 2 a 
(A 3 B 3 )/2L 2 , 



(ATI) 
(A.12) 
(A.13) 
(A.14) 
(A.15) 
(A.16) 



If we had considered a further cubic power of x in the ex- 
pansion (A.4), the previous expressions for the coefficients 
would have included another correcting term containing 
cubic power and cross-products of Ai and Bi multiplied 
by 1/3Lq, and the number of Fourier components would 
have been nine. 

We will assume that the absolute values of the co- 
efficients A 2 , B 2 are much smaller than those of A\ , B\ 
and A 3 , B 3 , that is, the second Fourier component is very 
small with respect to the first and third ones. We note 
that here we are not dealin g with the no nlinear oscilla- 
tor problem (e.g. Antonello 19944 |l994b|) . In the coeffi- 
cient of the second cosine and sine term, (A.6) and (A.12), 
the first elements, A 2 /L Q and B 2 /L a are, according to 
our assumption, small in comparison with the absolute 
value of the correcting terms which contain squares and 
cross-products of A x , A 3 , B X , B 3 . On the other hand, for 
the same reason the corrections of the coefficients of the 
first and third cosine and sine terms are small. In other 
words, while the first and third Fourier components are 
only slightly changed, we must expect a very different sec- 
ond component of the Fourier decomposed V light curve 
from that of the L light curve. This conclusion applies, of 
course, to any value of blending. 
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